Magnetic and thermal properties of 4f— 3d ladder-type molecular compounds 
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We report on the low-temperature magnetic susceptibilities and specific heats of the isostructural 
spin-ladder molecular complexes L2[M(opba)]3 ■ ^DMSO-yHbO, hereafter abbreviated with L2M3 
(where L = La, Gd, Tb, Dy, Ho and M = Cu, Zn). The results show that the Cu containing 
complexes (with the exception of La2Cu3) undergo long-range magnetic order at temperatures below 
2 K, and that for Gd2Cu3 this ordering is ferromagnetic, whereas for Tb2Cu3 and Dy2Cu3 it is 
probably antiferromagnetic. The susceptibilities and specific heats of Tb2Cu3 and Dy2Cu3 above 
Tq have been explained by means of a model taking into account nearest as well as next-nearest 
neighbor magnetic interactions. We show that the intraladder L-Cu interaction is the predominant 
one and that it is ferromagnetic for L = Gd, Tb and Dy. For the cases of Tb, Dy and Ho containing 
complexes, strong crystal field effects on the magnetic and thermal properties have to be taken into 
account. The magnetic coupling between the (ferromagnetic) ladders is found to be very weak and 
is probably of dipolar origin. 



PACS numbers: 75.30.-m; 75.30.Gw; 75.40.-s 
I. INTRODUCTION 



In today's search for smaller, faster, more selective and 
efficient products and processes, the engineering of well- 
defined spatial microarrangements of pure and composite 
materials is of vital importance for the creation of new 
magnetic devices. A possibility to assemble microstruc- 
tures in a controlled way is to use molecular-based ma- 
terials. The design of novel ferromagnetic molecular ma- 
terials is certainly among the stimulating subjectsi In 
such a context, molecular complexes based on transition 
metal ions are good candidates because the procedures 
for predicting the ferromagnetic nature of the interac- 
tion in this class of materials have become fairly well 
established^^ However, such an understanding is much 
less advanced when a lanthanide ion is involved. Mag- 
netic investigations concerning molecular materials con- 
taining lanthanide and transition metal ions have been 
overlooked until recently due to the weak interactions 



present and the complications added by the orbital con- 
tribution of lanthanide ions^ The case of Gd with a 
8 ^7/2 single-ion ground state and a quenched orbital con- 
tribution is special; most of the interest has been fo- 
cused on the Gd-Cu combination because of the very 
frequently found ferromagnetic character of the Gd-Cu 
interaction! 5 ! 6 ! 7 ! 8 ! 9 : 10 ! 11 : 12 : 13 : 14 Coupling those blocks in a 
3D lattice would produce molecular magnets and possi- 
bly ferromagnets. 

In an extended polynuclear complex, namely 
Gd 2 (ox)[Cu(pba)] 3 [Cu(H 2 0) 5 ] • 20H 2 O [hereafter 
abbreviated with Gd 2 Cri3(pba)] with pba = 1,3-propa- 
nediyl(oxamato), not only the ferromagnetic coupling 
of the Gd-Cu interaction was recently reported, but 
also the onset of a long-range magnetic orderii 5 ^ The 
crystal structure of this compound consists of layers of 
double-sheet polymers separated by water molecules. 
A puckered arrangement of Gd[Cu(pba)] units forms 
a two-dimensional honeycomb pattern connected by 
oxalato groups. Discrete [Cu(H 2 0)s] 2+ entities are 
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interspersed in the gap between the layers. From specific 
heat measurements, the onset of low-dimensional short- 
range order was found at temperatures around 1.5 K. By 
further lowering the temperature, the phase transition to 
three-dimensional ferromagnetic long-range ordering was 
observed at Tc — 1.05 K. This compound represents the 
first molecular-based ferromagnet involving lanthanide 
ions. 

Although long-range ferromagnetic order has thus been 
found in Gd2Cu3(pba), this material is isotropic due to 
the fact that both the Gd ions and the Cu ions are 
magnetically isotropic. A successful method to increase 
the intrinsic anisotropy is to synthesize intermetallic lan- 
thanide (L) and transition metal (M) compounds, where 
the L has an orbital contribution, that introduces the 
anisotropy that the transition metal is lacking. Such a 
strategy has of course been applied in research on metallic 
ferromagnets, such as Nd2Fei4B. The strong orbit cou- 
pling within L, on one side, and the intense L-M ex- 
change coupling has the net effect of polarizing the two 
sublattices in a direction, thus creating a strong uniaxial 
magnet, albeit, in a restricted temperature region^ The 
same method in trying to increase the uniaxial anisotropy 
has been applied in a series of molecular compounds 
based on L and M ions, by choosing the L that induces 
such anisotropy (first condition) , while expecting that the 
L-M interaction remains ferromagnetic and of similar in- 
tensity to that in the Gd2Cu3(pba) compound (second 
condition). The L anisotropy of the ground state de- 
pends on the crystal field acting on the orbital moment. 

In what follows, we shall describe the magnetic and 
thermal properties of a series of isostructural complexes 
based on lanthanide and transition metal ions. The gen- 
eral chemical formula of the compounds here studied is 
L 2 [M(opba)] 3 -a;DMSO-yH20 (hereafter abbreviated with 
L 2 M 3 ) where L = La, Gd, Tb, Dy, Ho and M = Cu, Zn, 
and opba stands for or</io-phenylenebis(oxamato), while 
DMSO stands for dimcthylsulfoxidc. In view of previ- 
ous work (Refii&) proving that the L-Cu exchange for 
light L substitutions is antiferromagnetic, thus not ful- 
filling the second necessary condition above mentioned, 
we have chosen to explore only the heavy L substitu- 
tions. Out of these, the most promising are the Tb and 
Dy substitutions, where there are already evidences of 
ferromagnetic coupling with Cu from susceptibility mea- 
surements, while no hint of such was found for Er, Tm or 
Yb>i& Recently, we have also reported, by means of spe- 
cific heat experiments, the onset of long-range magneti- 
cally ordered states in the Tb and Dy substituted com- 
pounds^ In this article we combine the thermal prop- 
erties of L 2 M 3 together with low-temperature suscepti- 
bility measurements, and explain the obtained results by 
means of theoretical calculations. 

Another interesting property of this series of isostruc- 
tural complexes is that they exhibit a spin-ladder struc- 
ture. Spin-ladders are low-dimensional magnetic quan- 
tum systems that consist of two or more strongly mag- 
netically coupled chains of spins and are thus intermedi- 



ate between one and two dimensional magnetic systems. 
The magnetic properties of such systems have received a 
renewed experimental and theoretical interest since the 
discovery of high-Tc superconductivity in ladder struc- 
tures provided by some cuprates like SrCu203f2& or in 
vanadyl pyrophosphate— or in LaCu02.5»^ The research 
on novel molecular complexes with ladder geometries is 
a very important task for a better understanding of the 
physics behind these complex systems^ 



II. STRUCTURAL AND EXPERIMENTAL 
DETAILS 

Information on the structure of the L2M3 series of com- 
pounds has been provided in detail in ReisM^. Here, 
we summarize the most important information in view of 
the analysis of the magnetic properties given below. 

The compounds crystallize very poorly and, conse- 
quently, a complete X-ray diffraction analysis could not 
be performed to determine the structure. Instead, the 
wide-angle X-ray scattering (WAXS) technique has been 
used to obtain structural information. Even though 
the WAXS technique cannot give the exact structure 
as compared to a full X-ray diffraction study, it has 
been successfully applied in several cases, for instance 
to structurally characterize inorganic polymers— i Fur- 
thermore, a complete structure determination by X-ray 
diffraction has been carried out on a related and very 
similar spin-ladder compound, namely, Tni2[Cu(opba)]3 ■ 
xDMF-yH 2 (a; « 10 ; y w 4). This crystal structure 
has therefore been utilized as a starting point to interpret 
the WAXS data for the L 2 M 3 compounds. 

The analyses strongly suggested that the structure of 
L2M3 is likewise based on a discrete, infinite ladder-like 
motif, as shown in Fig. 1. The sidepieces of the lad- 
der consist of L 2 [M(opba)]3 units developing along the 
b direction with an alternation of L(III) and M(II) ions. 
Rungs provided by the precursor [M(opba)] 2 ~, located 
between two L ions, connect the sidepieces. The L-M 
distance across an oxamato bridge is « 5 A, while the 
distance between two L ions across a M(opba) group 
is w 11 A. The shortest distance between two L ions 
(w 10 A) involves non-connected L ions belonging to two 
sidepieces of neighboring ladders. 

Each L ion is surrounded by seven oxygen ions, six 
of which belong to three opba ligands and one to a 
water molecule. The L coordination polyhedrons may 
be described as capped trigonal antiprisms. All the 
compounds are highly solvated with DMSO and H2O 
molecules, and the exact number of noncoordinated sol- 
vent molecules is not known accurately. Some of these 
molecules are easily removed. Therefore, there is some 
uncertainty in the molecular weight of the compounds 
and consequently in the absolute value of the molar sus- 
ceptibility and specific heat. 

The experiments done in the course of this work con- 
sisted primarily in measurements of three quantities: 
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magnetic susceptibility, magnetic moment and specific 
heat. Magnetic moment and susceptibility data down 
to 1.7 K were obtained with a commercial SQUID-based 
magnetometer with an ac-option. The ac-susceptibility 
of Gd2Cu3 in the temperature range 0.1K<T<3K was 
measured in Zaragoza with a mutual inductance tech- 
nique in a dilution refrigerator. The excitation ampli- 
tude was 10 mOe and the frequency / = 160 Hz. The 
signal was measured by means of a low-impedance ac- 
bridge, in which a SQUID was employed as a null de- 
tector^ The low-temperature susceptibilities down to 
5 mK of Tb2Cu3, Dy2Cu3 and H02CU3 were measured in 
Leiden with an ac-susceptometer directly mounted inside 
the mixing chamber of a dilution refrigerator. The fre- 
quency of the experiments was / = 520 Hz. The data in 
the overlap region between 1.70 K and 3 K were used to 
convert the low-temperature data from arbitrary into ab- 
solute units. The specific heat measurements on Gd20u.3, 
Gd 2 Zn 3 and Dy 2 Cu 3 in the range 0.2 K < T < 5.7 K were 
performed in Zaragoza by using an adiabatic calorime- 
ter cooled by adiabatic demagnetization, using the heat- 
pulse technique and Ge thermometry. The error on the 
specific heat has been estimated to be less than b%£& 
The specific heats of Tb2Cu3, H02CU3 and Ho2Zn 3 were 
measured down to 0.1 K in Leiden using a thermal relax- 
ation technique. The calorimeter was mounted in a dilu- 
tion refrigerator and connected to a cold sink through a 
calibrated heat link. All data were collected on powdered 
samples of the compounds. 



III. MAGNETIC PROPERTIES ABOVE 2 K 

Magnetic properties of systems containing lanthanide 
ions, such as L2M3, are significantly influenced by the in- 
teraction between a lanthanide ion and the surrounding 
ions in its direct environment. As a result, when a free 
lanthanide ion is placed in a crystal, its 2J+1 fold degen- 
eracy is partially lifted through electrostatic interaction 
between its /-electrons and the charges of the surround- 
ing ions. The multiplct is split into a number of states, 
which can appropriately be termed the crystal field (CF) 
states. Magnetic susceptibility measurements may offer 
valuable informations on the CF splittings. In this sec- 
tion we shall first present ac-susceptibility data above 
2 K of La 2 Cu 3 , Gd 2 Cu 3 and Gd 2 Zn 3 for which the CF 
effects can be neglected. Thereafter we shall briefly re- 
view and re-analyze in terms of the CF formalism the 
dc-susceptibility measurements on L2M3, with L = Tb, 
Dy, Ho and M = Cu, Zn, which already appeared in 
RefA 

The ac-susceptibility of Gd2Cu3 was measured with 
a frequency of 90 Hz and an amplitude of 4 Oe of the 
exciting field. The dc-measurements of Refi^ were per- 
formed in the temperature range 2-300 K, with a mag- 
netic field of 10 3 Oe. Diamagnetic corrections of the con- 
stituent atoms were estimated from Pascal's constants as 
—291 x 10~ 6 emu/mol. Due to the above mentioned un- 



certainty in the molecular weights, it was assumed that 
for each compound the maximum expected value for the 
Curie constant was reached at 300 K, and the experimen- 
tal data were rescaled accordingly when necessary. 

Hereafter we use the notation J, J' and J" for the 
exchange constants of the L-Cu, Cu-Cu and L-L inter- 
actions, respectively (negative values stand for antiferro- 
magnetic interactions) . 



A. La2Cua, GCI2CU3 and Gd2Zri3 

The magnetism of La2Cu3, where the lanthanide ion is 
diamagnetic, allows an estimate of the underlying Cu-Cu 
interaction. The single-ion magnetic properties of Cu(II) 
are fairly straightforward. Spin-orbit coupling causes the 
g values of the lowest Kramers doublet (S = 1/2) to 
lie in the range 2.0 to 2.3. Figure 2 shows the data of 
La2Cu3, which can be fitted to a Curie- Weiss law with 
8 = —0.2 K. This suggests that, even though the dis- 
tances in the ladder between next-nearest copper ions 
are as long as « 10 A, they still interact with each other 
with a weak antiferromagnetic superexchange coupling. 
Taking S — 1/2, 6 — -0.2 K and z = 4 for the number of 
nearest neighbors for each copper ion, the mean-field ex- 
pression for the Curie- Weiss temperature^ provides the 
value of J' /k,R = —0.2 K for this Cu-Cu interaction. 

Gadolinium(III) has a 8 S 7 / 2 ground state; the or- 
bital contribution is almost entirely quenched and very 
isotropic g values close to the free electron value are 
found. For temperatures above 2 K, the in-phase ac- 
susceptibility of Gd2Cu3 follows a Curie- Weiss law with 
= 2.3 K and C — 16.4 emu K mol -1 , indicating the fer- 
romagnetic nature of the compound (Fig. 2). In the limit 
T 3> 6, the paramagnetic susceptibility can be described 
by 



C 
T 



2 C*Gd + 3 Ccu 
f : 



(3.1) 



where Ccd = ^Me//(Gd)/3fcB is the Curie constant 
of the two gadolinium ions and Cq u = N /(Cu) /3&B 
is the Curie constant of the three copper ions. Taking 
S Gd = 7/2, 5 Cu = 1/2 and g Gd = g Cu = 2.00, one 
obtains 



c = 2 a 



Gd 



3 a 



Cu 



16.8 emu K mol 



which is in satisfactory agreement with the experimen- 
tal value (16.4 emu K mol -1 ). The same mean-field anal- 
ysis used above gives the value of J//cb = 0.5 K for the 
Gd-Cu exchange coupling. The origin of the ferromag- 
netic nature of the Gd-Cu interaction is further discussed 
in Section IVAl 

When the copper in Gd2Cu3 is substituted by diamag- 
netic Zn, then the susceptibility of Gd2Zn3 could be ex- 
pected to follow the Curie law for two isolated Gd ions. 
Nevertheless, Fig. 2 shows that the Curie behavior (with 
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the expected value of C — 15.7 emu K mol -1 ) is no longer 
maintained below 10 K, since a slight decrease of the \F 
product from 15.7 down to 14.8 emu K mol^ 1 at 2 K is 
observed. The data are better described by a Curie- Weiss 
law with 6 = —0.1 K. Such a deviation may be due to 
a weak, antiferromagnetic Gd-Gd interaction of dipolar 
origin. Since the Gd moment is large, dipole-dipole inter- 
actions could show up in the susceptibility at tempera- 
tures as high as 10 K, given the nearest-neighbor Gd-Gd 
distance of 10 A. Of course it may also be attributed to 
a weak antiferromagnetic superexchange coupling acting 
between the Gd ions. The mean- field analysis provides 
the value of J"/^b = —6 x 10~ 3 K for this interaction, 
on basis of the obtained #-value. 



B. Tb 2 Zn 3 , Dy 2 Zn 3 and Ho 2 Zn 3 

As mentioned, for a lanthanide ion with a non-zero or- 
bital moment, such as Tb, Dy or Ho, the effect of the 
crystal field on the magnetic levels has to be taken into 
account. In order to study the CF splittings for L2M3, we 
consider the susceptibility measurements for the M = Zn 
compounds to avoid (or reduce) the complication of mag- 
netic exchange interactions. To simplify the analysis we 
consider the lanthanide ions in a cubic geometry, so that 
the number of crystal field parameters is minimized. As 
a further simplification, we consider an octahedral co- 
ordination for each lanthanide ion instead of the seven 
coordination found for Tni2[Cu(opba)]3 • xDMF-yH^O 
(x ~ 10 ; y ~ 4), which crystallizes in orthorhombic 
geometryiSiS^ Although these simplifications may seem 
somewhat insatisfactory, the lack of a precise knowledge 
of the actual coordination parameters (distances, angles) 
for the various L ions in L2M3 does actually not justify 
more sophisticated treatments. Within this approxima- 
tion, the prediction of the type and the relative order of 
the CF split levels is greatly simplified. We use the valu- 
able paper of Lea, Leask and Wolf 30 (hereafter referred 
to as LLW) which predicts both order and type of level in 
cubic symmetries for each value of the angular momen- 
tum and for all ratios of the crystal field parameters. 

The crystal field acting on the lanthanide ion in octa- 
hedral geometry can be represented by the Hamiltonian 

Hcf = B A (Ol + 50 4 ) + B 6 (O 6 - 210%), (3.2) 

where the O™ terms are angular momentum operators. 
The parameter B n is related to the strength of the crystal 
field components 

B n = a n A° n (r n ), 

where a n is the operator equivalent factor, a, (3, 7 for 
n = 2,4,6, respectively. are parameters relating the 
nth degree potential on the ion from the ionic charges 
of the lattice, and (r n ) are the expectation values for 
the distances of the 4/ electrons. The product A°(r n ) is 



determined to the first order by geometric factors and in 
principle can be calculated. 

LLW have determined the eigenfunction and eigen- 
value solutions for the Hamiltonian equation for appli- 
cable values of the angular momentum and the full range 
of Bi/B^ ratios. They have constructed diagrams for 
the eigenvalues of the crystal field levels of each angular 
momentum in terms of a quantity x, defined as 

x F{A)B A _ F(4)/3Ag(r 4 ) _ b 4 
1 - \x\ ~ F(6)B 6 ~ F(6) 7 A°{r 6 ) _ V 

where — 1 < x < 1, and F(4) and F(6) are multiplica- 
tive factors. Then, the susceptibility can be easily calcu- 
lated for a given angular momentum and a gyromagnetic 
ratio. In the remainder of this section, we discuss the 
susceptibility data^ of L2Zn3 (with L = Tb, Dy, Ho) 
compounds and we explain them in terms of the LLW 
scheme. Only the temperature region below 40 K will 
be considered, since only the contribution of the lowest 
lying levels needs to be taken into account. These ap- 
proximations are allowed in view of the large size of the 
level splittings compared to the low-temperature region 
of interest for our experiments (mainly below 10 K). This 
is confirmed by the specific heat analyses given below, 
where we shall show that only few CF levels are actually 
involved in the magnetic ordering processes. 

Let us start with Tb 2 Zn 3 . Compounds of non-Kramers 
lanthanide ions, such as Tb(III), often have a singlet elec- 
tronic ground state separated by an energy A from an 
excited singlet state. In a cubic field, the 7 F 6 state of Tb 
is split into two singlets (ri and I^), one non-magnetic 
doublet (r 3 ) and three triplets [r 4 , and r| 2) ]. 

The measured susceptibility of Tb2Zn3 is shown in 
Fig. 3. The data at low temperatures are well described 
by two singlets split by A ~ 0.2 K. In the LLW scheme 
for octahedral coordination, such a situation corresponds 
to |x| ~ 0.45 where the ground state can be either Ti or 
T2 (see RefiSi) . The assignment will be confirmed by the 
specific heat analysis presented below. For this simple 
case of two separated singlets, to improve the quality of 
the fit, we used ultimately a model in which a weak an- 
tiferromagnetic exchange interaction acting between the 
terbium ions has been taking into account. The fit pro- 
vides the value of J" /k-g = —1 K for this interaction. 
The model used is explained in detail in the Appendix 
and employed in Section IV Bl to explain the very-low- 
temperature magnetic and thermal properties of Tb2Cu3. 

For the case of dysprosium(III), the crystal field splits 
the 6 i?i5/2 ground state into two Kramers doublets (r 6 

and Tf) and three quartets fig, and T^) (see 

Refi22>). In octahedral coordination with |x| » 0.45, as 
estimated for the previous case of Tb, one anticipates a 
doublet as the ground state (Tg or T7) with the higher 
state being the other doublet (F7 or r 6 , respectively). 
The fit to the Dy 2 Zn 3 data (Fig. 3) is in agreement with 
a low-lying doublet being either Tg or T7 separated by 
A = 13.7 K from or Tq, respectively. Also this as- 
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signment will be confirmed by the analysis of the specific 
heat (see below). 

Holmium(III) is a non-Kramers ion, with a 5 Is ground 
state. The crystal field splits the ground state into one 
singlet (Ti), two non-magnetic doublets [r^ and T^] 
and four triplets [T^ , rf\ if } and if ]. The LLW di- 
agram (see Refi^S) is complicated, with many level cross- 
ings, which make simple choices difficult. For octahedral 
coordination with \x\ w 0.45, as in the Tb and Dy cases, 

(2) (2) 

one expects a T y 3 ' ground state with T\ ' the first ex- 
cited state within a distance of not more than ~ 20 K. 
The fit to the Ho2Zn3 data (Fig. 3) is in agreement with 
this, yielding the value of A = 8.3 K for the T^-T^ 
separation. 

To summarize, the susceptibility data presented here 
for the Tb, Dy and Ho containing compounds, as well as 
the specific heat analyses presented below, can be rea- 
sonably accounted for by assuming cubic symmetry and 
octahedral coordination for the lanthanide ions. The con- 
sistency of the analyses is derived from the fact that the 
factor x could be taken at approximately the same value 
(sa 0.45) for each compound. This is what one would 
expect for a series of isostructural compounds in which 
the lanthanide ions most probably have the same local 
coordination. 



IV. VERY LOW-TEMPERATURE PROPERTIES 

In the following sections, we present the thermal prop- 
erties of the L2M3 compounds in the temperature range 
0.1 to 10 K, and their magnetic properties down to 5 mK. 
We first analyze the lattice contributions to the measured 
heat capacities, and then describe the results inferred 
from susceptibility, magnetization and magnetic specific 
heat for each compound of the series. 



A. Phonons contribution 

Figure 4 shows the collected zero-field specific heat 
data of the L2M3 compounds as a function of tem- 
perature. As a first step, we have estimated the lat- 
tice contribution for each compound by fitting the high- 
temperature data to a sum of a Debye lattice contribu- 
tion (Ci/R = (3 T 3 ) and a high-temperature limiting tail 
(oc T~ 2 ) of the magnetic anomaly. The results for the j3 
parameter are given in Table I. 

The WAXS analysis has shown that the L2M3 com- 
pounds have similar structures (see Section [H] and ref- 
erences therein). Consequently, one expects also a sim- 
ilar lattice contribution for each of them. Indeed, the 
estimated (5 values for five of the compounds studied 
fall into a fairly narrow range between 1.1 xlO -2 and 
3.2 xl0~ 2 K~ 3 . Since the coefficient (3 depends on the 
third power of the Debye temperature, such a limited 
variation appears quite acceptable. The only exception 



is Dy 2 Cu3 for which f3 — 8.8 x 10~ 2 K~ 3 , a value much 
higher than the others. The reason is presently not 
known. It is important to note that, however, the mag- 
netic specific heat for Dy2Cu3 obtained after subtracting 
the phonon specific heat, appears to yield a magnetic en- 
tropy that agrees with the value expected for a doublet 
ground state. 

The magnetic contributions to the specific heats of all 
the L 2 M 3 compounds discussed in the following sections 
were obtained by subtracting the above mentioned lattice 
contributions. 

B. Gd 2 Cu 3 and Gd 2 Zn 3 

We have already seen ( Section IIII Afl that for temper- 
atures above 2 K, the in-phase susceptibility of Gd2Cri3 
follows a Curie-Weiss law with 9 — 2.3 K, indicating 
the ferromagnetic nature of the compound. The low- 
temperature in-phase susceptibility of Gd2Cu3 is de- 
picted in Fig. 5 as a function of the temperature. A 
sharp peak is observed at Tc = 1-78 K which is proba- 
bly due to a transition to long-range magnetic order. A 
large value of the susceptibility is found at the maximum 
(65.9 emu/mol). This behavior is typical for a powdered 
sample of an isotropic ferromagnetic material in which 
demagnetization effects become important. The theoret- 
ical value estimated for the susceptibility of a ferromag- 
netic sample of Gd2Cu3 at Tc is XextC^c) = N^ 1 « 
(70 ±20) emu/mol, where N is the demagnetizing factor- 
calculated taking into account an ellipsoidal approxima- 
tion for the geometry of the sample. Within the error 
this is equal to the experimental result. 

An interesting feature, that is discussed later (Sec- 
tion IV A|) , is the in-phase susceptibility variation with 
the temperature for T < Tc (Fig. 5). Just below Tc, 
the susceptibility goes through a minimum and further 
lowering the temperature reveals a rounded anomaly for 
0.1K<T<1K. 

Figure 6 shows the variation of the experimental mag- 
netization M versus the field H for the temperature 
T = 1.70 K, i.e. just below the critical temperature 
(Tc = 1.78 K). No hysteresis effect is observed. Addi- 
tional evidence for the ferromagnetic coupling is gained 
by comparison of the magnetization measurements to the 
calculated behavior for the uncoupled case. If the mag- 
netic centers were all uncoupled, M would vary according 
to 

M = 2Ng Gd /iBSGdB 7 /2(riGd) + SNgcuHBScuBi^ivCu) 

where 5 7 / 2 (?7Gd) and B 1 / 2 (r]Cu) are the Brillouin func- 
tions for Gd and Cu ions, respectively. The saturation 
magnetization would be given by 

M = {2g Gd S Gd + 3g C uScu)Nfi B = 17N^ B - 

The calculated magnetization for T = 1.70 K as de- 
rived from Eq. (|4.1|l is displayed in Fig. 6 together with 
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the experimental curve. It can be observed that the ex- 
perimental data saturate at the predicted value and they 
lie above the expression of Eq. j|4.1|l in the whole tem- 
perature range. Such a behavior clearly indicates that 
the predominant interaction, i.e. the Gd-Cu interaction 
through the oxamato bridge, is ferromagnetic. 

Let us now discuss the magnetic specific heat data 
by presenting first the result for the Gd2Zn3 compound 
(Fig. 7). A sharp increase appears when lowering the 
temperature below 1 K. Due to incomplete achievement 
of the ordering process in the temperature range of our 
experiments, an analysis of the experimental entropy is 
not possible. 

Figure 8 shows the magnetic molar specific heat C m /R 
versus the temperature T for Gd2Cu3. Two clear fea- 
tures can be observed: a distinct A-peak at Tc = (1.78 ± 
0.02) K confirming the assignment of a long-range or- 
dering process, and a rounded maximum around 0.8 K. 
We note that the rounded anomaly occurs in the same 
temperature range as that observed in the susceptibil- 
ity experiments (Fig. 5). Figure 8 shows also that the 
magnetic contribution at T > Tc is quite large and ex- 
tends up to high temperatures, indicating the presence 
of important short-range ordering effects most probably 
related to the low dimensionality (no contributions from 
CF splittings are expected for Gd ions). 

In order to calculate the entropy, we use the relation 

/>oo 

S/R= / (C m (T)/R-T) dT, (4.2) 
Jo 

together with the experimental values of C m in the 
temperature range 0.2 K < T < 5.7K. For the ex- 
trapolation down to K we assume a 3D ferromagnetic 
spin-wave type contribution (oc T 3 / 2 ) and on the high- 
temperature side a T~ 2 dependence for C m . The calcu- 
lated entropy gives a value of S/R « 9 In 2. This corre- 
sponds to the maximum expected value (= 3 In 2 + 2 In 8) 
evidencing that both the Gd and Cu magnetic ions par- 
ticipate in the ordering process. As already said, the 
relatively large values of the specific heat above Tc may 
be associated with low-dimensional fluctuations within 
the ladders. This can be clearly seen by plotting the en- 
tropy variation as a function of the temperature (inset of 
Fig. 8). The calculated entropy variation above Tc gives 
(Soo — Sc) ~ 2.7 R In 2, which is a rather large fraction in 
comparison with theoretical values of three-dimensional 
models of ferromagnets^ All this agrees with the low- 
dimensional ladder-type magnetic structure of these ma- 
terials. 



C. Tb 2 Cu 3 

The low-temperature susceptibility of Tb2Cu3 is de- 
picted in Fig. 9. The abrupt change of the in-phase sus- 
ceptibility at 1 K is ascribed to a transition to a magnet- 
ically ordered state, which is obvious also from the spe- 



cific heat curve shown in Fig. 10. The maximum value is 
93.6 emu/mol at 1.0 K. This value is almost of the same 
order as expected for a ferromagnetic material in which 
demagnetization effects become important. In fact, an 
estimate of the demagnetizing factor gives the value of 
« 280 emu/mol for the susceptibility of an isotropic fer- 
romagnetic Tb2Cu3 sample at the maximum, assuming 
a cylindrical approximation of its shape. At the low tem- 
perature side of the anomaly, the susceptibility decreases 
sharply, reaching a value of 5.0 emu/mol at 10 mK. The 
specific heat data are plotted in Fig. 10 and indicate a 
transition temperature of Tc — (0.81 ± 0.01) K. In the 
susceptibility data we see that fluctuations in the order- 
ing process show up as a peak in the out-of-phase sus- 
ceptibility x" centered at 0.7 K (inset of Fig. 9), whereas 
the in-phase susceptibility x' shows an inflection point 
(maximum in dx' /dT) at that same temperature. The 
x" is in arbitrary units because it is null above 0.8 K 
and, thus, it was not possible to scale to data measured 
above 1.7 K with the commercial magnetometer. The ob- 
served behavior is not uncommon for a low-dimensional 
system which undergoes a phase transition to long-range 
magnetic order, for instance to an antiferromagnetic ar- 
rangement of the ferromagnetic ladders, at Tc = 0.81 K 
as deduced from the inflection point of x' ■ The relatively 
large value of the susceptibility at 1 K may suggest the 
presence of a dominating ferromagnetic interaction, prob- 
ably associated with the intraladder Tb-Cu interaction. 
We will return to this point in Section IV Bl Moreover, 
the very low value reached by the in-phase susceptibil- 
ity in the limit of very low temperature, clearly indicates 
the large anisotropy of the terbium ions in Tb2Cu3, as 
already pointed out in Section fill Bl for Tb2Zn3. The 
anisotropy may also be responsible for the fact that the 
measured value of the x' a t Tc is lower than the expected 
limit for an isotropic ferromagnetic sample. Since a pow- 
der is measured, the value is a directional average. 

The main feature of the magnetic specific heat of 
Tb 2 Cu 3 , depicted in Fig. 10, is the A-peak at T c = 0.81 K 
indicating the transition to a long-range ordered state. 
The bump observed above Tc and centered around 1.5 K 
is probably associated with low-dimensional fluctuations 
within the ladders, since a contribution from excited CF- 
levels is not expected in this temperature range (see Sec- 
tion ^^Bj). The data at the low-temperature side show 
a slight upward curvature. Below 0.4 K, the data obey 
the law C m T 2 /R = 0.05 K 2 . This contribution is proba- 
bly coming from the hyperfine splitting of the magnetic 
levels of the Tb nuclei. Subtracting the calculated hyper- 
fine specific heat from the measured specific heat data, 
we estimate the remaining entropy. In order to calcu- 
late the magnetic entropy of Tb2Cu3, we carry out the 
integration in Eq. 14. 2|) using the experimental values of 
C m /R in the temperature range 0.1 K < T < 5 K. After 
extrapolation down to K with an exponential function, 
the calculated entropy gives a value of S/R ss 5 In 2. Tak- 
ing into account that Cu has spin 1/2, which corresponds 
to S/R = In 2 per atom, the experimental entropy con- 
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tent indicates that Tb has an effective spin 1/2, corre- 
sponding to a lowest lying doublet in the temperature 
region of the magnetic ordering process. In Section UlI Bl 
we have shown that for Tb2Zn3 the lowest energy lev- 
els are two singlets with a separation of A = 0.2 K. For 
Tb2Cu3 we expect a similar crystal environment and, as 
a consequence, a similar energy scheme. Nevertheless, 
the separation A is relatively small in comparison with 
the thermal energy in the region of the magnetic ordering 
process. We thus conclude that we can apply the approx- 
imation of effective spin 1/2 for the Tb ion. In agreement 
with the pronounced high-temperature specific heat tail 
observed in Fig. 10, only about 15% of the entropy is left 
below Tc- This is another indication of the large amount 
of short-range magnetic order that is probably associated 
with the low-dimensionality of the ladder. 



D. Dy 2 Cu 3 

Figure 11 shows the low-temperature ac-susceptibility 
of Dy 2 Cu3. As in the case of Tb2Cu3, and as evi- 
denced from the specific heat experiment presented in 
Fig. 12, clear evidence of a transition to a magnetically 
ordered phase is found. The in-phase susceptibility has 
a sharp peak centered at 0.8 K with a maximum value of 
116.4 emu/mol. Below the peak, it decreases very sharply 
and remains almost constant at 10.0 emu/mol down to 
10 mK. As for Tb2Cu3, the maximum of the experimen- 
tal susceptibility is not far below the expected theoretical 
limiting value (~ 280 emu/mol) for an isotropic ferro- 
magnetic sample. The out-of-phase susceptibility also 
shows a sharp peak centered at 0.77 K and it is zero be- 
low 0.4 K and above 1 K (inset of Fig. 11). At 0.77 K an 
inflection point occurs in the temperature dependence of 
the x' ■ The transition temperature deduced below from 
the specific heat is Tc = (0.75 ± 0.01) K. Similar as for 
Tb2Cii3, the data point toward a type of antiferromag- 
netic long-range ordering, with very small value for the 
antiferromagnetic coupling ( Jaf) between the ferromag- 
netic ladders. Since for an antiferromagnet the x & t Tc is 
inversely proportional to Jaf, when the latter becomes 
small enough, x a t Tc can reach the demagnetizing limit. 
It is then difficult to distinguish between fcrro- or antifer- 
romagnetic coupling between the ferromagnetic chains. 
The peak in the x" found in both compounds may be 
due to a weak-ferromagnetic moment which arises when 
the antiferromagnetic ordering is accompanied by some 
degree of spin canting. Again, as for Tb2Cu3, the very 
low value of the in-phase susceptibility for T — > indi- 
cates large anisotropy of the Dy ions. 

Figure 12 shows the magnetic specific heat, plotted 
as C m /R versus T, for the Dy2Cu3 compound, where 
C m /R is the molar specific heat and T is the tempera- 
ture. The prominent spike below 1 K is identified with 
the A anomaly indicating the onset of a phase transition 
to a long-range ordered state. The peak is very sharp 
and allows an accurate determination of the critical tem- 



perature as T c = (0.75 ± 0.01) K. The analysis of the 
magnetic entropy content shows that only a very small 
portion (about 15%) is obtained below Tc and that, as 
for Tb 2 Cu3, the lanthanide has an effective spin 1/2 in 
the ground state, once more agreeing with the crystal 
field analysis ( Section |lll B|) . 



E. H02C113 and Ho2Zri3 

The low-temperature susceptibility of H02CU3 is shown 
in Fig. 13. Also for this compound, evidence of an an- 
tiferromagnetic long-range ordering is observed. The in- 
phase susceptibility has a maximum of 27.5 emu/mol at 
0.12 K (much lower than for Tb2Cu3 and Dy 2 Cu3). It 
sharply decreases down to 8.9 emu/mol at 30 mK, and 
remains nearly constant by further lowering the tempera- 
ture down to 5 mK. The out-of-phase susceptibility shows 
a peak around 40-100 mK, in which temperature range 
also an inflection point, centered at (50 ± 10) mK, is seen 
in the in-phase susceptibility. 

The molar specific heat of Ho 2 Cu3 is depicted in Fig. 14 
together with that of Ho2Zn3 for comparison. The spe- 
cific heats of H02CU3 and Ho 2 Zn3 are seen to overlap 
for T > 0.8 K, showing a broad rounded anomaly with 
a maximum at 3.5 K (Fig. 14). This is probably due 
to crystal field splitting of the holmium levels. In Sec- 
tion ^^B] we have seen that the susceptibility of Ho2Zn3 

(2) 

down to 2 K may be explained in terms of a doublet (T 3 ) 

ground state separated by 8.3 K from a triplet (T^) ex- 
cited level. Taking into account a further splitting of 
these degenerate levels, the Schottky contribution aris- 
ing from this configuration can be easily calculated. A 

(2) 

nice agreement with the data is obtained for T 3 further 

split by 2 K, and separated by 11 K from T\ , further 
split by 5 K each level (Fig. 14). 

For H02Z113, a sharp increase of the specific heat by 
lowering the temperature is observed below 0.3 K. This 
contribution is apparently coming from the hyperfine 
splitting of the magnetic levels of the Ho nuclei. In this 
temperature region, by fitting the data to C m j R = aT~ b , 
we obtain a — 0.14 and b — 1.3, and not the expected 
value of b = 2, for a magnetic anomaly high-temperature 
tail. The reason for that is probably related to prob- 
lems of thermal contact between spin system and the 
lattice, since we found the specific heat at low tempera- 
ture to become increasingly dependent on the measuring 
time we used in the experiments. In fact, both specific 
heats of H02CU3 and Ho 2 Zn3 have been measured with 
the relaxation technique. If the measuring time is not 
long enough for the sample to achieve thermal equilib- 
rium, between the spin system and the phonon bath, part 
of the electronic spin contribution will simply be miss- 
ing. Also for Ho 2 Cu3 an upward curvature is observed 
by lowering the temperature below 0.7 K. By compari- 
son with the Ho2Zn3 results, we see that in the specific 
heat of Ho 2 Cu3 there is an extra contribution besides 
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the hyperfine splitting of the Ho nuclei. This may read- 
ily be explained by a weak Ho-Cu coupling that, together 
with the dipole-dipole interladder interaction, shows up 
at very low temperature. Consequently, only the onset 
of a phase transition is observed. The low-temperature 
limit of our setup did not allow to determine the ordering 
temperature that has to be below 0.1 K (judging from the 
X-data). In combination with the susceptibility results 
(Fig. 13), we conclude that the ordering temperature Tc 
for Ho 2 Cu 3 is between 0.04 and 0.10 K. 

Due to the incomplete achievement of the magnetic 
ordering processes at the lowest temperature, the analysis 
of the entropy contents for H02CU3 and Ho2Zn3 is not 
possible. 



V. COMPARISON WITH LADDER MODELS 
A. Gd 2 Zn 3 and Gd 2 Cu 3 

For Gd2Zii3, the sharp increase of the specific heat be- 
low 1 K (Fig. 7) may reflect the onset of a 3D ordering 
process, which takes place at a temperature lower than 
obtained in the experiment. Even though the Gd 2 Zn 3 
compound is more poorly crystallized than the Cu con- 
taining compounds^ its measured specific heat clearly 
underlines the relevance of the Gd-Cu interaction, when 
compared with the results found for Gd2Cu3. The ab- 
sence of the copper at M sites reduces or inhibits the or- 
dering temperature of the Gd sublattice below our lowest 
experimental temperature. It is clear that the Cu ion in 
Gd 2 Cu3 is much more effective than the Zn ion in trans- 
mitting magnetic exchange interaction along the Gd- 
oxamato-Cu-oxamato-Gd superexchange pathway that 
gives rise to the net ferromagnetic coupling. 

In order to explain the paramagnetic susceptibility of 
Gd 2 Cu 3 , let us consider the model proposed by Georges 
et al^ for ladder-type double chains which has been 
successfully applied to describe the paramagnetic sus- 
ceptibility of the Gd 2 (ox)[Cu(pba)] 3 [Cu(H 2 0) 5 ] • 20H 2 O 
compound. This model, which is based on the standard 
transfer matrix method, takes into account the presence 
of copper quantum spins and gadolinium classical spins, 
in a ladder-like arrangement similar to the L 2 M3 com- 
pounds, where each gadolinium ion interacts isotropically 
with two neighboring copper ions. The approximation is 
the classical treatment of the gadolinium spin, which is 
allowed because of the high spin value of the ion. 

With the aid of this model, we have analyzed the exper- 
imental thermal dependence of the susceptibility in the 
paramagnetic regime. The same Lande factor g = 2.00, 
determined in the high-temperature limit, has been at- 
tributed to all cations. A unique coupling constant J 
has been introduced for all the Gd-Cu isotropic interac- 
tions. The best fitting of the experimental susceptibility 
results is obtained in the range 2.5K<T<20K. Fig- 
ure 5 shows the calculated curve with J/k-Q = 0.74 K. 
As expected, the positive J value refers to a ferromag- 



netic exchange coupling, in agreement with the positive 
paramagnetic Curie temperature = 2.3 K. It is worth- 
while to note that the J value obtained here does not 
differ substantially from the estimate of 0.5 K given in 
Section fill Al deduced by simple mean field analysis. At 
temperatures below 2.5 K, the three-dimensional order- 
ing mechanism becomes apparent and, consequently, the 
experimental data can no longer be described with the 
above (paramagnetic) model. 

It is well established that an array of isotropic spins 
with dimensionality two or less will not present long- 
range ordering^ The experimental specific heat data 
have shown the presence of short-range order for tem- 
peratures far above Tc already, and we have associated 
it with extended magnetic correlations along the spin lad- 
ders. The observed long-range ordering at Tc = 1.78 K 
implies the existence of interactions between adjacent 
ladders. The structure of Gd 2 Cu3 does not present any 
pathway between adjacent ladders that involve atomic 
contacts suitable for the propagation of magnetic ex- 
change. Consequently, it seems reasonable to assume 
that the driving force for the magnetic long-range order 
is indeed the dipolar interaction. 

Another interesting feature in Gd 2 Cu3 is the maximum 
observed at about 0.8 K, and thus below Tc — 1.78 K, 
in both the susceptibility (Fig. 5) and the specific heat 
(Fig. 8). To explain such a feature, we may consider the 
presence of two types of interchain coupled systems in 
Gd 2 Cu3, e.g. with slightly different packing of the lad- 
ders within the crystal structure. Then one polytype has 
Tc = 1-78 K and the other has a lower critical tempera- 
ture of 0.8 K. Unfortunately, due to the lack of a precise 
knowledge of the structure and to the fact that the ex- 
periments have been performed on powdered samples, a 
detailed study of the ordering process could not be car- 
ried out. 

For the sake of completeness we wish to mention briefly 
the nature of the mechanism of the Gd-Cu interaction. 
In preceding works^ 1 ^ the ferromagnetism of this inter- 
action was attributed to the coupling between the 4/-3d 
ground configuration and the excited configuration aris- 
ing from the metal-metal charge transfer configurations. 
The latter is associated with the 3d — > 5c? process: an 
electron is transferred from the singly-occupied orbital 
centered on copper toward an empty orbital centered on 
gadolinium. In such a mechanism, J is given by 

J = Y / ^3 d S/^U 2 -~S 2 )] l (5.1) 

z=l 

where ftsd-sd is the transfer integral of the 3d —> 5d 
process, S is the energy gap between S = 3 and 5 = 4 
excited states arising from the 4/ 7 5d 1 electron-transfer 
configuration and U is the energy cost of such a transfer. 
The summation applies to the five bd gadolinium orbitals. 

It is interesting to compare, qualitatively, the here ob- 
tained coupling constants J of Gd 2 Cu3 {J/ks — 0.74 K) 
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and of Gd 2 (ox)[Cu(pba)]3[Cu(H 2 0) 5 ] • 20H 2 O, which is 
a long-range ordered ferromagnet for T < 1.05 K±£ii& 
and has J/k B = 0.40 The difference may in fact 
be understood on the basis of the above mechanism. 
In the latter compound, the Gd ions adopt a distorted 
capped square antiprism coordination of oxygen atoms; 
six out of a total of nine oxygen atoms are provided by 
the three bridging Cu(pba) groups, two by the oxalato 
ligand, and the last one by a water molecule^ In the 
former compound, due to the absence of the oxalato lig- 
and, the 5d orbitals of the gadolinium ions are distributed 
in a less distorted environment. This fact, together with 
slightly shorter Gd-Cu distances for Gd 2 Cu3 (see Refs^ 
and2^ for comparison), gives higher values of f3$d-zd for 
Gd 2 Cu3. Consequently, according to Eq. (|5.1|) . the cou- 
pling constant J should indeed be larger for Gd 2 Cu3 than 
for Gd 2 (ox)[Cu(pba)] 3 [Cu(H 2 0) 5 ] • 20H 2 O, as found. 



B. Tb 2 Cu 3 , Dy 2 Cu 3 and Ho 2 Cu 3 

Our measurements show Tb 2 Cu3, Dy 2 Cu3 and 
Ho 2 Cu3 to be strongly anisotropic systems with long- 
range order phase transitions below 1 K. As for Gd 2 Cu3, 
the low-dimensional structure of the compounds implies 
that the driving force for the long-range magnetic or- 
der is the dipolar interaction between the ferromagnetic 
chains. However, for Tb 2 Cu3 and Dy 2 Cu3, evidence for 
an even larger amount of short-range order, associated 
with fluctuations within the ladder structure, is found 
above Tc- From the specific heat data, we have deduced 
that Tb 2 Cu3 and Dy 2 Cu3 may be described in terms of 
an effective spin 1/2 at low temperatures. We note that 
an effective spin 1/2, arising from a low- lying doublet, 
is often found for lanthanidc compounds of dysprosium 
and terbium with low magnetic ordering temperature^ 
At higher temperatures the behavior is strongly modified 
by the presence of low-lying excited states. 

To explain the magnetic and thermal properties of 
Tb 2 Cri3 and Dy 2 Cri3, a model has been developed, which 
is based on the standard transfer matrix technique, using 
the Ising approximation and assuming spin 1/2 for the L 
and the Cu ions. The details of the model can be found in 
the Appendix. As revealed from the susceptibility mea- 
surements, we have to expect a ferromagnetic coupling 
J for the L-Cu interaction. Next-nearest neighbor in- 
teractions are also taken into consideration. In fact, we 
have seen in Section fill Al that a weak antiferromagnetic 
interaction is acting between the copper ions in La 2 Cu3. 
Consequently, we may assume a similar Cu-Cu interac- 
tion to be present in L 2 Cu3 for L = Tb and Dy. We shall 
denote by J' the Cu-Cu interaction. Similarly, we have 
also taken into account a weak antiferromagnetic cou- 
pling J" acting between the L ions, as suggested by the 
susceptibility data for Gd 2 Zn3 (Section IIII A|) . Accord- 
ing to our notation, a positive coupling constant implies 
a ferromagnetic coupling. The presence of next-nearest 
neighbor interactions in ID molecular systems based on 



lanthanide ions is not new and several examples can be 
found in the literature^ For the case of Tb 2 Cu3, the 
crystal field energy separation A of the two lowest sin- 
glets of the non-Kramers Tb ion has also been included 
in the calculation. 

To improve the reliability of the fits for Tb 2 Cu3 and 
Dy 2 Cu3, we independently use the model in La 2 Cu3 to 
get the estimate of J', and in Tb 2 Zii3 and Dy 2 Zn3 to 
get the estimates of J". We then fit the susceptibilities 
and specific heats of Tb 2 Cu3 and Dy 2 Cu3 with J' and 
J" fixed at the values obtained precedently with only 
the coupling J and the effective gyromagnetic constant 
g as free parameters. The fits to the susceptibility and 
specific heat of Tb 2 Cu3 are shown in Figs. 9 and 10, re- 
spectively. The fits provide the values of J/k B = 5.6 K, 
J'/k B = -1.0 K, J"/k B = -1.0 K, A = 0.2 K and 
an effective gyromagnetic constant of g — 18.0 for the 
Tb ion. Similarly, the results for the susceptibility and 
specific heat of Dy 2 Cu3 are shown in Figs. 11 and 12, 
respectively. In this case, the fits provide the values 
of J/k B = 4.6 K, J'/ka = -I -0 K, J"/k B = -1.0 K, 
and effective g = 19.6. It can be seen that the exper- 
imental behaviors above the 3D ordering temperatures 
are satisfactory reproduced. Only the specific heat of 
Tb 2 Cu3 is not so well accounted for (Fig. 10). Proba- 
bly, this comes from the influence of low-lying excited 
levels, making the Hamiltonian of this system less Ising- 
like than for Dy 2 Cu3. Moreover, the large g values found 
for both systems imply large magnetic moments which 
are moreover coupled ferromagnetically into chains, so 
that dipolar interactions may be strong enough to con- 
tribute also above the 3D ordering temperatures. Con- 
sequently, the model here presented may fail to give a 
detailed explanation in this temperature region. Nev- 
ertheless, the following conclusions can be drawn. The 
compounds Tb 2 Cu3 and Dy 2 Cu3 are very similar in be- 
havior and are strongly anisotropic. We note that large g 
values of order 20, found for both Tb 2 Cu3 and Dy 2 Cu3, 
have been also reported for other Dy and Tb based com- 
pounds^ For both cases, the L-Cu interaction is pre- 
dominant and ferromagnetic. This interaction is slightly 
stronger for Tb 2 Cu3 than for Dy 2 Cu3, in line with the 
analysis reported in Ref. 1 — . Consequently, also the 3D or- 
dering temperature is higher for Tb 2 Cu3 (Tc = 0.81 K) 
than for Dy 2 Cu3 (Tc = 0.75 K). For both compounds, 
the presence of next-nearest neighbor interactions have 
to be taken into account. These interactions are much 
weaker than the L-Cu interaction and are antiferromag- 
netic in nature. In conclusion, the results here obtained 
suggest that Tb 2 Cu3 and Dy 2 Cu3 order ferromagneti- 
cally within the ladders. This ordering process is accom- 
panied by weak dipolar and probably antiferromagnetic 
interladder interactions, that, together with strong crys- 
tal field effects, lower the susceptibility below Tc- 

For Ho 2 Cu3, we have already reported that the Ho- 
Cu interaction has to be very weak. Indeed, crystal 
field effects and dipolar interactions are mainly responsi- 
ble for the magnetic properties at low temperature (Sec- 
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tion IIV E|) . Consequently, also the 3D ordering temper- 
ature is the lowest one in comparison with Tb 2 Cu3 and 
Dy 2 Cu 3 . 

The results obtained for the magnetic interaction pa- 
rameters are summarized in Table II for the Cu contain- 
ing compounds and in Table III for the Zn containing 
compounds. 

VI. CONCLUDING REMARKS 

In the previous sections we have focused attention on 
the physical properties of spin-ladder molecular magnets 
containing lanthanide and transition metal ions. To- 
gether with the originality of their crystal structure, the 
L2M3 compounds present interesting magnetic features 
such as the onset of long-range orderings for L = Gd, 
Tb, Dy, Ho and M = Cu. To the best of our knowledge, 
Gd2Cu3 is the lanthanide and transition metal ions based 
ferromagnet with the highest long-range ordering critical 
temperature Tc = (1.78 ±0.02) K so far reported. More- 
over, as estimated from specific heat and susceptibility 
measurements, Tb2Cu3 has a magnetic ordering temper- 
ature of T c = (0.81 ± 0.01) K, while Dy 2 Cu 3 orders at 
T c = (0.75 ±0.01) K, and Ho 2 Cu 3 has a T c between 0.04 
and 0.10 K. These molecular based magnets are the first 
ones obtained with lanthanide ions other than gadolin- 
ium. 

The very pronounced quasi-one-dimensionality of the 
magnetic structure implies that the 3D ordering is driven 
by the dipolar interaction acting between ladders to- 
gether with the intraladder superexchange interaction. 
We have also reported that the intraladder L-Cu inter- 
action is the dominant one and it is ferromagnetic for 
L = Gd, Tb and Dy. The influence on the magnetic 
properties of the weaker antiferromagnetic next-nearest 
neighbor L-L and Cu-Cu interactions have also been an- 
alyzed. A remarkable point is the key role of the copper 
ions in these complex systems. They transmit the mag- 
netic exchange interaction between the L ions across the 



ladder. As a proof of this, we have seen that, if Cu is sub- 
stituted by the non-magnetic Zn, the L-L interaction is 
inhibited and the critical temperature of the 3D ordering 
process is strongly diminished. Finally, the effects on the 
susceptibility and specific heat of the crystal field split- 
tings of the magnetic energy levels of the lanthanide ions 
in L2M3 have been discussed and rationalized, assuming 
the same symmetry and coordination for the lanthanide 
ion in each compound of the series. 



APPENDIX: ISING MODEL FOR QUANTUM 
SPIN LADDER 

The model here reported has been developed to explain 
the magnetic and thermal properties of the Tb 2 Cu3 and 
Dy 2 Cu3 compounds. 

As a first step, we assume only a coupling between the 
next-nearest copper (Cu) and lanthanide (L = Tb, Dy) 
ions which is ferromagnetic (in the following notations, 
this means J > 0). 

For Tb and Dy, we assume an effective low- 
temperature spin 1/2 and Ising coupling due to the crys- 
tal field anisotropy of the lanthanide ion. For Cu, we 
assume spin 1/2 as well. The effective Hamiltonian is 
then: 



H = -J ^2(<TiSj) - {^gaWOi +^gs^SAH z , 

<=1 J=l 

(6.1) 

where cr, and Sj are the projections along the z axis of 
the spins of the copper and lanthanide ions, respectively 
(Fig. 15). The Lande factor of the copper ions is assumed 
to be g a — 2, while the one of the lanthanide ions is a 
fitting parameter. The Zeeman terms in Eq. (|6.1J) extend 
over all N a copper and N$ lanthanide ions. 

We build the transfer matrix on the unit cell 



<7i(Si + S[) + a 2 (S 2 + S! 2 )+ i <7 3 (Si + S' 2 ) + \ o±{S x + S' 2 ) 



T(S 1 ,S 2 :,S[,S' 2 ) = ex P Y J 

x exp ji (3gs^H z (Si + S[ + S 2 + S 2 ) + f3g a ^i B H z ^cti + <t 2 + i <t 3 + ^ 04 



where j3 — l/k^T. The partition is expressed in terms 
of the 4x4 transfer matrix 



Z(T,H Z ) = Y,T(S 1 ,S2;S{,S' 2 )T(S' 1 ,S 2 ;S?,S 2 '). 



...T(s[ N \si N) ;S 1 ,S 2 )=Tr(T N ) 
Therefore, by defining 

/ (T, H z ) = lira — In Z = In A 

where A is the largest eigenvalue of T, we obtain numer- 
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ically the magnetic susceptibility and the molar specific 
heat by 

According to the susceptibility data of La2Cu3 (Fig. 2), 
we assume a slight antiferromagnetic next-nearest neigh- 
bor coupling J' between the copper ions. Similarly, ac- 



cording to the susceptibility data of Gd2Zn3 (Fig. 2), we 
expect an antiferromagnetic next-nearest neighbor cou- 
pling J" between the lanthanidc ions. In both cases, the 
data are compatible with couplings not exceeding 1 K. 
By adding these second neighbors couplings we get the 
complete scheme depicted in Figure 15. 

The transfer matrix is now 



T(Si,a 3 ,S 2 ;S' 1 ,a 4 ,S' 2 ) = exp < (3 J 

<Tl,(T2 ^ 



a^St + S[) + <j 2 (S 2 + S' 2 )+ X - a 3 (S 1 + S' 2 ) + X - a 4 (Si + S' 2 ) 



1 



SiS[ + S 2 S' 2 + — {S\S 2 + S[S 2 ) 



x exp |/3J'(ai + cr 2 )(a 3 + cr 4 ) + (3J" 
x exp | i Pgs^HziSi + S[ + S 2 + S' 2 ) + (3g a ^ B H z ^ + u 2 + 1 a 3 + ^ u A 



with J > 0, J' < and J" < 0. This is now an 8 x 8 
matrix but the process is the same. 

For the terbium case, that is a non-Kramers ion, we 
allow the spin states to be non degenerate by adding to 
the Hamiltonian a term E(Si) such that E(+^) = — 
and E(— |) = — ^> where A is the gap. The best fits for 
Tb 2 Cu3 are displayed in Figs. 9, 10 and for Dy 2 Cu 3 in 
Figs. 11, 12. 
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Gd2Cu3 


Tb 2 Cu 3 


Dy2Cu3 


H02CU3 


Gd 2 Zn 3 


Ho 2 Zn3 





3.2 


1.5 


8.8 


1.2 


1.1 


2.9 



TABLE I: Experimental lattice contributions (Ci/R = (3 T 3 ) 
as estimated from the data in Fig. 4. The values f3 are ex- 
pressed in (xl0~ 2 K~ 3 ). 





La 2 Cu 3 


Gd 2 Cu 3 


Tb 2 Cu 3 


Dy 2 Cu 3 


Ho 2 Cu 3 


T C (K) 




1.78(2) 


0.81(1) 


0.75(1) 


0.04^0.10 


0(K) 


-0.2 


2.3 








J/k B (K) 




0.74 


5.6 


4.6 




J'/ka (K) 


-0.2 




-1.0 


-1.0 




J"/k B (K) 






-1.0 


-1.0 




A(K) 






0.2 (s-s) 




« 11 (d-t) 






2 


18.0 


19.6 




Model 


MF 


Georges 


App. 


App. 


Schottky-LLW 



TABLE II: Experimental results obtained for the Cu contain- 
ing compounds. Tc is the long-range magnetic ordering tem- 
perature; 8 is obtained from the Curie- Weiss law (Section ing : 
J, J' and J" are the exchange constants for the L-Cu, Cu-Cu 
and L-L interactions, respectively (negative values stand for 
antiferromagnetic interactions); A is the separation between 
the ground state and the first excited state (included also is 
the type of state: s=singlet; d=doublet; t=triplet); g^ is the 
Lande factor for the L ion; the models used are MF (mean 
field) or LLW ( Section lHTBl or Georges f Section lV^l i or App. 
(Appendix and Section FVBI or Schottky fSection llV El . 





Gd 2 Zn 3 


Tb 2 Zn 3 


Dy 2 Zn 3 Ho 2 Zn 3 


9(K) 


-0.1 






J"/kB (K) 


-6 x 10" 3 


-1.0 




A (K) 




0.2 (s-s) 


13.7 (d-d) 8.3 (d-t) 


Model 


MF 


LLW- App. 


LLW LLW 



TABLE III: Experimental results obtained for the Zn contain- 
ing compounds. The same notation of Table II is used. 
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LIST OF FIGURES 

1. Relative dispositions of the ladders: (a) view of two 
neighboring ladders in the be plane; (b) schematic 
view of the ladders projected in the ac plane. The 
lanthanide ions are located on the sides of the lad- 
der and occupy each corner, while each transition 
metal ion is between two lanthanide ions along the 
sides and in the rungs. 

2. TOP: Experimental molar susceptibility of La2Cu3. 
In the inset: the same data plotted as \F vs T . 
BOTTOM: Inverse of the in-phase susceptibilities 
of Gd 2 Zn 3 (•) and Gd 2 Cu 3 (o). In the inset: yT 
vs T for Gd 2 Zn3. 

3. TOP: Experimental inverse susceptibility for 
Tb 2 Zn3. The solid line is the calculated inverse 
susceptibility taking into account A = 0.2 K and 
J"/k B = -1 K for the Tb-Tb interaction. For 
explanations, see text. CENTER: Measured and 
calculated inverse susceptibility for Dy 2 Zii3. BOT- 
TOM: Measured and calculated inverse susceptibil- 
ity for Ho 2 Zn 3 . 

4. Experimental specific heats of L 2 M3. 

5. Experimental molar in-phase susceptibility of 
Gd 2 Cu3 versus temperature and theoretical esti- 
mation for J/k-Q — 0.74 K (see Section IV Al for 
explanations). 

6. Field dependence of the magnetization and hystere- 
sis loop at T = 1.70 K for Gd 2 Cu 3 . The dotted 
line represents the Brillouin functions as calculated 
from Eq. I|4.1|l . 

7. Magnetic specific heat vs temperature for Gd 2 Zn3. 
Solid curve is a guide to the eye. 

8. Magnetic molar specific heat vs temperature for 
Gd 2 Cu3. Curve a is the ferromagnetic spin wave 
contribution (oc T 3 / 2 ); curve b is the high tempera- 
ture limit (oc T~ 2 ) of the magnetic anomaly tail. In 
the inset: temperature dependence of the magnetic 
entropy for Gd 2 Cu3. 

9. Experimental low-temperature in-phase suscepti- 
bility vs temperature for Tb 2 Cu3 together with the 
calculated susceptibility (for explanations see Sec- 
tion ^^J. In the inset: Experimental out-of-phase 
susceptibility. 

10. Magnetic molar specific heat of Tb 2 Cu3. The dot- 
ted line is the calculated hyperfine contribution of 
the Tb ions. The solid line is the calculated specific 
heat due to the low-dimensionality of the ladder 
(for explanations see Section IV Bp . 

11. Experimental low-temperature in-phase suscepti- 
bility vs temperature for Dy 2 Cu3 together with the 



calculated susceptibility (for explanations see Sec- 
t ion IV B|) . In the inset: Experimental out-of-phase 
susceptibility. 

12. Magnetic molar specific heat of Dy 2 Cu3. The solid 
line is the calculated specific heat due to the low- 
dimensionality of the ladder (for explanations see 
Section EU. 

13. Experimental molar in-phase and out-of-phase (in- 
set) susceptibilities vs temperature for Ho 2 Cu3. 
Solid curves are guides to the eye. 

14. Magnetic molar specific heat of Ho 2 Cu3 (•) and 
Ho 2 Zn3 (o). The solid line is the calculated Schot- 

(2) 

tky contribution due to the splitting of the r 3 ' and 

(2) 

T\ ' levels of the holmium ions as shown in the Fig- 
ure. In the inset the same data in log-log scale. 

15. Ladder structure considered by the model discussed 
in the Appendix. 
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Fig. 3, M. Evangelisti et al., PRB 
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Fig. 9, M. Evangelisti et al., PRB 
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Fig. 10, M. Evangelisti et al., PRB 
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Fig. 11, M. Evangelisti et al, PRB 
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